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We use light-cone QCD sum rules to calculate the pion-photon transition form factor, 
taking into account radiative corrections up to the ncxt-to-next-to-leading order of per- 
turbation theory. We compare the obtained predictions with all available experimental 
data from the CELLO, CLEO, and the BaBar Collaborations. We point out that the 
BaBar data are incompatible with the convolution scheme of QCD, on which our predic- 
tions are based, and can possibly be explained only with a violation of the factorization 
theorem. We pull together recent theoretical results and comment on their significance. 
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1. Introduction 

For many years now the production of a pion by the fusion of two photons has 
attracted the attention of theorists and experimentalists. Theoretically, the proc ess 
l*{li)l*{l2) ^ T^'^ip) can be treated within the convolution scheme of QCeUJ^ by 
virtue of the factorization theorem which allows one to treat the photon-parton 
interactions within perturbative QCD. while all binding effects are separated out 
and absorbed into the pion distribution amplitude (DA). This latter ingredient has 
a nonperturbative origin and can, therefore, not be computed within perturbative 
QCD. One has to apply some nonperturbative approach to derive it or reconstruct 
it from the data. A widespread framework to calculate static and dynamical non- 
perturbative quantities of hadrons is provided by QCD sum rules with locaf^l or 
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nonlocal condcnsatesPl This latter method was employed by us in collaboration 
with A. P. Bakulev to derive a pion DA that is able to give good agree- 

ment with various sets of data pertaining to various pion observables, e.g., the 
electromagnetic form factoJ^, the pion-photon transition form factoJ3, diffractive 
di-jets productionjSl etc. 

While the process with two off-shell photons is theoretically the most prefer- 
able, experimentally, another kinematic situation is more accessible, notably, when 

one of the photons becomes real, as probed by the 

CELLCP and the CLEO 



CollaborationJl^. Such a process demands more sophisticated techniques in or- 
der to take properly into account the hadronic content of the real photon. Indeed, 
first Kh odja mirian,!^ then Schmedding and Yakovlev,22I used light-cone sum rules 
(LCSR)fP^ to analyze the CLEO data, a method also applied by BMS up to the 
next-to-leading-order (NLO) level of QCD perturbation theory.L^^ 
Especially the high-precision CLEO datcPl 

on F'''^ give the possibility to 
verify the pion DAs quantitativeljff^H^H^I^. It was found that the best agreement 
with the CLEO data is provided by pion DAs which have suppressed endpoints 
a; = 0, 1, like those belonging to the "bunch" determined by BMS irPwith the help of 
QCD sum rules with nonlocal condensates. Note that the endpoint suppression is a 
sui generis feature of the nonlocality of the qu ark condensate and is controlled by the 
vacuum quark virtuality « 0.4 GeV^ .^^El All these approaches attempt to reverse 
engineer the pion DA from its (first few) moments. In the BMS approactP the 
first ten moments have been calculated, from which the corresponding Gegenbauer 
coefficients a„ with n = 0, 2, . . . , 10 were determined. It turns out that all coefficients 
with n > 4 are negligible, so that the proposed model DA has only two coefficients: 
a2 and 04. 

Recently, wjl^ extended this type of calculation to the NNLO of QCD pertur- 
bation theory by taking into account those radiative corrections at this order which 
are proportional to the /3o-function. More specifically, we used the hard-scattering 
amplitude of this order, computed before in Ref. I17( in order to determine the spec- 
tral density within the LCSR approach mentioned above. In addition, we refined the 
phenomenological part of the SR by using a Breit-Wigner ansatz to model the me- 
son resonances. Below, we report about the main results of this analysis and further 
discuss what conclusions can be drawn by comparing the obtained predictions with 
all the available experimental data. We focus attention on the new BaBar data,^^^ 
which turn out to be incompatible with our predictions, indicating a violation of 
collinear factorization in QCD. We perform a detailed comparison of these data 
with the theoretical expectations and some proposed scenarios to explain them. 
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2. Pion-photon transition form factor F''' in QCD 

The transition form factor F'' describes the process 7*((?i)7*((?2) """(p) and 
is given by the following matrix element {—Qi = > 0, — = 9^ > 0) 



d^xe-^^^-^TT^ip) I T{j,{z)MO)} I 0) 



te^,^^qUl-F^''''^Q\q^). (1) 



Provided the photon momenta are sufficiently large Q^,q^ ^ (where the hadron 
scale is set by the p- meson mass Wp), the pion binding effects can be absorbed into 
a universal pion distribution amplitude of twist-two. Then, one obtains the form 
factor in the form of a convolution by virtue of the coUinear factorization.UJ^ 



'{Q\q^) = ^(Q^ q^^4■,x) ® (pI'H^; mI) + o (q-^) 



(2) 



Here the pion DA ip"^' represents a parametrization of the pion matrix element 
at the (low) factorization scale /ip, whereas the amplitude T, describing the hard 
parton subprocesses, can be calculated in QCD perturbation theory: T = To + 
agTi -\- a^^T2 ■ ■ where as = Q!s/(47r), and with O {Q~^) denoting the twist-four 
contribution. In leading order (LO) of the strong coupling and taking into account 
the twist-four term explicitly, one 



Nf 



dx 



ip^^\x;^j.l) 
Q^x 



q^x 



dx 



[Q'^x -t- q^x)'^ 



with Nf = ^/tt and x=l 



- X 



The pion DA of twist-two. Lp\ \ is defined by 



q{z)-f,^^C{z,Q)qm<P)) „ = *Pp/. / dxe"(^-^Vi'H^,Ml), 



(3) 



(4) 



where C{z,Q) = Vexp{ig A^{T)dT'^) is a path-ordered exponential to ensure 
gauge invariancc. The second term in Eq. ^ represents the twist-four contribution, 
which is becoming important for small and intermediate values of Q^. The pion 
DA of twist-four. (p^^\ is an effective onJISl in the sense that it is composed from 
different pion DAs of twist-four. In our analysis it is taken in its asymptotic form. 
The parameter is det erm ined from the matrix element {'K{P)\gsdGa^'y'^u\0) ~ 
iS'^Up^,, and is estimateciHl to be (5^(1 GeV^) = 0.19±0.02GeV2. Estimates for the 
twist-four term, based on the rcnormalon approach,^^ have been considered in the 
last entry of Ref . IHl 

On the other hand, also the evolution of tp^n\x,^p) with /i| is controlled by 
a perturbatively calcula ble e volution kernel V, following the Efremov-Radyushkin- 



Brodsky-Lepagc (ERBLpJ^ equation 



fi^^^ix; = [as (x, y) + al V^^\x, y) + ..)® ^(y; ^J?) (5) 



d^jL 

T/(°) ® V'n = 2Cf v{n) ■ ^. 
v{n) = l/{n -|- l){n 



M^)^6xx C(3/2)(.x-s); (6) 



l/2 + 2(V;(2)-V^(n + 2)); tl;{z) ^ -HT{z)) . (7) 



4 S. V. Mikhailov, N. G. Stefanis 



Here, {ipnix)} are the Gcgenbaucr harmonics, which constitute the LO eigenfunc- 
tions of the ERBL equation, v(n) being the corresponding eigenvalues. Then, one 
has 

^i^Hx; fi^) = M^) + an(M') V'n(a;), (8) 

n=2A,... 

where the coefficients {a„} evolve (in LO) with /i^ and have specific values for each 
pion DA model (a compilation of the coefficients of various proposed models can be 
found in Refs. [Til [71) . 

The radiative corrections to the hard amplitudes in NLO, encapsulated in Ti, 
have been computed in Ref. [201 More recently, the /3-part of the NNLO amplitude 
T2, i.e., Pq ■ T2_/3, was also calculated. It is instructive to discuss the structure of 
this result at the scale /Xp = especially in view of further considerations in Sec. 
131 



, (9) 



where L{y) = In [(Q^y + g^y)/^'!] ■ The first term, CpTg^^^ is the /3o-part of the 
NNLO coefficient function and represents, from the computational point of view, 
the most cumbersome element of the calculation in Ref.[T7l The next term originates 
from the NLO coefficient function T^^-* using one-loop evolution of a,. The third 
term appears as the /3o-part of the two-loop ERBL evolution (see Eq. while 
the last term, which is proportional to V^^'^ , stems from the combined effect of the 
ERBL-evolution and the one- loop evolution of as- 



3. Light Cone Sum Rules for the process 7*(Q^)7(q^ ~ 0) ^ 7r° 

The transition form factor, when one of the pho tons bec omes quasi real (g^ — *■ 0), 
has been measured by different Collaborations El 1 01 1 ^ 1 However, this kinematics 
requires the modification of the standard factorization formula Eq. ((J) in order to 
take into account the long-distance interaction, i.e., the hadronic content of the 
on-shell photon. A viable way to reach this goal is to employ the method of LCSRs, 
which are based on a dispersion relation for F'^ in the variable g^, viz., 

i^^*-^*^ (Q2, q') = r ■ (10) 

Jo s + q 

The key element in this equation is the s pect ral density p{Q^,s) ~ 
iHi [^7*7*^^(^2 „g)j fQj. ^jjj^jj jjjj^jjg ^ijg ansatffH p = pP'^Q^ s)6'(so - s) + 
p^'^ {Q^ , s)9{s — So), where the "physical" (ph) spectral density p^^ serves to ac- 
commodate the hadronic content of the photon (below an effective threshold sq) by 
means of the transition form factors F'^ ^'^ of vector mesons, notably, p or w. 



pP\Q^ s) = y2/pF^*^-(02) . S{s ~ ml) 



(11) 
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On the other hand, p^"^ embodies the partonic part of the LCSR and derives from 
Eq. ^ via the relation p^'^{Q'^,s) = ^ [F^' '''''' (Q'^ , -s)]. The pion-photon tran- 
sition form factor F'' ''"^((5^,0) can be expressed in terms of p^"^ having recourse 
to quark-hadron duality in the vector channel: 



Im (F^'^'^^iQ^-s)) 



ds 



So 

1 Im 

^ Jo 



ds , 



(12) 



with So — 1-5 GeV'^ and A/^ denoting the Borel parameter in the interval (0.5 — 
0.9) GeV^. The LO spectral density p^''^ has been determined in Ref. [11] using for 
pi'i'^ Eq. dSl). 

Partial results for the NLO spectral density p^^^ in the leading-twist approx- 
imation have been given in Ref. [121 while the general solution pn\Q'^,s) = 
Jsi [(Ti ® Vn) (Q^ -s)] was obtained in Ref. [Uto read 



Pn^ (x; ^1) 



(13) 



with 



-3 [1 - v'^in)] + ^ - In^ (^) + 2v{n) 

6 \x/ \X /ip 



-Cf 2 ^ [Gni + w(n) • 6„i] {x) 
1=0,2,... 



V'n(2;) 
(14) 



where w° (n) = l/(n+l)(n + 2) — 1/2 and Gni , bni are calculable triangular matrices 
(see Ref. [in] for details). Note that the spectral density pi^"* allows one to obtain 
F^QSYiiQ^) ^'^^ number of the Gegenbauer harmonics in the expansion of (pir- 

Employing this approach, predictions at the NLO for Q^F^Q^^iQ^) were 
derived jlSl using a variety of pion DAs: the asymptotic (Asy) one, the CZ model,l21 
and the BMS "bunch"^ (see left panel of Fig. [1]). It was found that the radiative 
corrections arc important, being negative and contributing up to -17% at low and 
moderate values. Let us recall the main results of this CLEO-data analysis, refer- 
ring to Refs.[T4][7lfor a full-fledged discussion and more details. The CLEO datcP^ 
were processed in terms of a error ellipses and the results are displayed in the right 
panel of Fig. [1] around the best fit point (+). In this graphics, the BMS "bunch" of 
pion is shown as a slanted green rectangle, while the vertical dashed and solid 

lines denote the estimates for 02 (related to the second moment of (Pt^) of two recent 
lattice simulations in Refs. [211 ES] respectively. Several other models, explained in 
Refs. [14] [7] are also shown. The main upshot is that wide pion DAs, like the CZ 
model, are excluded at the 4cr level, whereas the asymptotic DA, and others close 
to it, are also excluded at the level of at least 3cr. As one sees from this figure. 
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only cndpoint-supprcsscd pion DAs of the BMS type are within the la error elhpse 
of the CLEO datEp^ and simuhaneously in agreement with the mentioned lattice 
constraints. 



Q'F.,^,{Q'-) [GcV] 





Fig. 1. Left: Predictions for Q^^Lcgp(Q^) using the CZ model (upper dashed red line), the BMS 
"bunch" (shaded green strip), and the Asy DA (low dashed black line) in comparison with the 
CELLO (diamonds) and the CLEO (triangles) data. Right: CLEO-data constraints on ^lcsr('3^) 
in the ( ap . a d) plane at the scale /^^ = (2.4 GeV)^ in terms of error regions around the BMS best-fit 
point using the following designations: la (thick solid green line); 2a (solid blue line); 3a 

(dashed-dotted red line) . Two recent lattice simulation s!^ ^ are denoted, respectively, by vertical 
dashed and solid lines together with predictions of QCD sum rules with nonlocal condensates 
(slanted green rectangle), * — BMS model,l51* — asymptotic DA, ■ — CZ DA.El 



The inclusion of the main, i.e., /3o-proportional NNLO contribution, in ^"lcsr 
proceeds via the dispersion integral in (fT^ .P^ The technical problem is how to 
obtain the contributions to pt^^/^) from terms with various powers of the logarithms 
L{y) in Eq. The outcome of this calculation turns out to be negative, like 
the NLO contribution, and about -7% (taken together with the effect of a more 
realistic Breit-Wigner (BW) ansatz for the meson resonances in Eq. ([TT]) ) at small 

~ 2 GeV^. The size of this suppression decreases rather fast to -2.5% with 
increasing > 6 GeV^. The net result is a slight suppression of the prediction for 
the scaled form factor (see Fig. [5]). 



4. Confronting NNLO LCSR results with the BaBar data 

In the preceding section we have shown in detail that the CLEO data are strictly 
incompatible with wide pion DAs and demand that the endpoints a; = 0, 1 are 
stronger suppressed than in the asymptotic DA. Surprisingly, the new data of the 
BaBar Collaboratioiff^ on the pion-photon transition form factor are in contra- 
diction with this behavior. More specifically, these data, which extend from inter- 
mediate up to high momenta in the range A < < AO GeV^, show a significant 
growth with for values above ^ 10 GeV^. Indeed, the corresponding data points 
lie above the asymptotic QCD prediction v^/t and continue to grow with up 
to the highest measured momentum. This behavior of the BaBar data is clearly 
in conflict with the coUinear factorization. This in turn means that, as we argued 
in Ref. [HI the inclusion of the NNLO radiative corrections cannot reconcile the 
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Fig. 2. Left: Predictions for Q^^lcsr ^it^ the NNLO/3 corrections (together with the BW 
resonance model) included (solid lines) and without them (dashed lines). They correspond to 
selected pion DAs: CZ model — upper (red) lineP BMS-model — middle (green) linepl and Asy 
DA — lower (black) line. The experimental data are as in Fig. [T] Right: The shaded (green) 
strip denotes the form-factor predictions derived for the EMS "bunch" .1^ The dashed (blue) line 
represents a dipole-form interpolation of the CLEO data. 



BaBar data with perturbativc QCD. This is true for any pion DA that vanishes 
at the endpoints a; = 0, 1 (see Fig. [3]). From Tabled] it becomes clear that also a 
wide pion DA, like the CZ model, though it lies above the asymptotic prediction 
shows exactly the same scaling behavior above ^ 10 GeV^ as all discussed pion 
DAs, which is not even surprising because it also vanishes at the endpoints x — ,1. 
The conclusion is that — contrary to the statements of the BaB ar CoUaboartiorM^ 
— also the CZ DA cannot reproduce all BaBar Indeed, from Fig. [3] 

one sees that in the region of the CLEO data, the CZ model fails also with respect 
to the BaBar data points at the 4(t level, whereas above ~ 20 GeV^, and up to the 
highest measured value of Q^, 40 GeV^, it fails again, because instead of growing 
with it scales. 



Table 1. xLf ^^'^ Asymptotic (Asy), BMS, and CZ DAs. 



Pion 
DA 


BaBar & 


z CLEO data BaBar all data 


BaBar only highest 10 data 

X-ndf 


Asy 


11.5 


19.2 


19.8 


BMS 


4.4 


7.8 


11.9 


CZ 


20.9 


36.0 


6.0 



In summary, 

(i) The combined effect of the negative NNLO^, i.e., suppressing, radiative correc- 
tions and the enhancing effect of using a Breit-Wigner model for the vector-meson 
resonances finally amounts to a moderate overall suppression of Fj^csr(Q^) 
the range of momentum transfer 10-40 GeV^ jlS'^ probed in the BaBar experiment. 

(ii) The growth with of the scaled form factor, measured by the BaBar Collab- 
oration, cannot be attributed to the hadronic content of the real photon because 
this is a twist-four contribution that is rapidly decreasing with increasing Q^. 

(iii) It is impossible to get enhancement of the form factor within the QCD collinear 
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Fig. 3. Predictions for F^q^{Q^) calculated with the following pion DAs: Asy — lower solid 
line . BM S "bunch" — shaded green strip, and the CZ model — upper solid red lin e. The BaBar 
dataisJ are shown as diamonds with error bars. The CELLclSl and the CLEcHSI data are also 
shown with the designations used in Fig. [T] The displayed theoretical results include the NNLO^ 
radiative corrections and the BW model for the meson resonances. The horizontal dashed line 
marks the asymptotic QCD prediction v^/vr. 



factorization using pion DA models which have a convergent projection onto the 
Gegenbauer harmonics and, hence, vanish at the endpoints x = 0, 1 (cf. the val- 
ues in Table [T|). It seems (see next section) that exactly the violation of this feature 
may provide an explanation of the BaBar data. 



5. BaBar data — heuristic explanations 

Despite the claims by the BaBar CoUaboratiorff^ that their data are in agreement 
with QCD, such an explanation is outside of reach at present. Hence, once is forced 
to look for alternative explanations. There have been several proposals to explain 
the anomalous behavior of the BaBar data, among others, e.g., Refs. EH [26l [27l 
[28l [29l [30l We restrict attention to one sort of such proposals which assumes that 
the pion DA may be "practically flat", hence violating the coUinear factorization 
and cntaihng a (logarithmic) growth of Q'^F'^'^'^ {Q"^) with Q^. One haPSl (p-Z = 
0.53 GeV2) 



V2U 



1 



dx . 



(15) 



Another optioiJ^ gives instead (m sa 0.65 GeV, cf. Eq. ^) 

q2^7*7x ^ V2/^ ip^ix, Q) 



-dx 



x + 



(16) 



with (p^(a;,/io) = (1 - N)&xx with N w 1.3 and /iq = 6-0 8 GeV. Equation 
P5p can be compared with the available experimental 

datPMm 

using for them 



Pion transition form factor at the two-loop level vis-a-vis experimental data 9 



the phenomcnological fit (A w 0.9 GcV, 6 ~ —1.4) 



A2 



A2 + Q' 



+ b 



A2 



A2 + Q2 



by means of a Xndf criterion, given in Tabled One observes — diagonal from the 



Table 2. X^df phenomcnological fit and the fit with a 

fiat pion DA, like Eq. JTSj — numbers in parentheses. The diagonal 
elements give the best-fit values of X^df ^^^'^^ fi^ the corresponding 
line parameters. 





CELLO& CLEO 


BaBar 


CELLOiSCLEO 


0.48 (1.22) 


7.8 (15.8) 


BaBar 


10.8 (3.5) 


1.8 (1.8) 



BaBar (CELLO&CLEO) entry to the CELLO&CLEO (BaBar) one — that it is not 
possible to fit simultaneously the CELLO/CLEO data and the BaBar data with the 
same accuracy using these parameterizations. 

This interpretation is given further impetus, when we consider the BaBar data as 
being two "experiments" BaBarl and BaBar2 — see Fig.[5l One sees from Table 3, in 
terms of Xndf ' that the flat-DA scenario cannot describe both BaBar "experiments" 
simultaneously with the same accuracy. 



. 25 
0.2 
. 15 



Q' [GeV^] 



Table 3. xL/ ^^'^ fiat-DA fit. 
Best-fit values on the diagonal. 




♦ 

BaBarl 


■ 

BaBar2 


BaBarl 


3.3 


0.33 


BaBar2 


3.5 


0.26 



5 10 15 20 25 30 35 



Fig. 4. Split BaBar data: ♦ - BaBarl, ■ - BaBar2. 



6. Conclusions 

We have studied in detail the pion-photon transition form factor using light-cone 
sum rules and including QCD radiative corrections up to the two loop level. We also 
took into account twist-four contributions. It has been our goal to derive predictions 
for Q^F''^ ^ using several pion distribution amplitudes that can be compared with 
the available experimental data. Our results have been deduced within the con- 
volution scheme of QCD for distribution amplitudes that vanish at the endpoints 
X = 0, 1. They turn out to be unable to match the data of the BaBar Collaboration 
for momenta beyond 10 GeV^, which grow with increasing . We analyzed this 
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behavior and argued that proposed scenarios, which make use of fiat pion distribu- 
tion amplitudes, cannot match the high-Q^ BaBar data and those of the CLEO and 
the CELLO Collaborations simultaneously, because the latter demand distribution 
amplitudes that vanish at the endpoints. 
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